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We present a comprehensive first-principles study of the ballistic transport properties of low di- 
mensional nanostructures such as linear chains of atoms (Al, C) and carbon nanotubes in presence of 
defects. A novel approach is introduced where quantum conductance is computed from the combi- 
nation of accurate plane-wave electronic structure calculations, the evaluation of the corresponding 
maximally-localized Wannier functions, and the calculation of transport properties by a real-space 
Green's function method based on the Landauer formalism. This approach is computationally very 
efficient, can be straightforwardly implemented as a post-processing step in a standard electronic- 
structure calculation, and allows to directly link the electronic transport properties of a device to 
the nature of the chemical bonds, providing insight onto the mechanisms that govern electron flow 
at the nanoscale. 

PACS numbers: 71.15.Ap,72.10.-d,73.63.-b 



I. INTRODUCTION 



The field of nanotechnology has undergone a remark- 
able growth in the last few years. This development 
has been fueled by the expectation that unusual prop- 
erties of matter^™* which become evident as the dimen- 
sions of the structural components of a device shrink 
under ~10-100 nm, may be exploited. Indeed, the 
challenges for future developments involve continuous 
shrinking of the physical dimensions of devices and 
attainment of higher speeds. The drive to produce 
smaller devices has led the current research towards 
new forms of electronics, in which nanoscale objects and 
molecular devices replace the transistors of today's sili- 
con technology4i£*£i£ Experiments have been performed 
to directly measure charge transport properties of hy- 
brid metal-(bio)molecular gystems^^iSiiSiiiiiSiiiii and 
carbon-based aggregates such as fullerenes and nan- 
otubes ,MiiSiiSiiLiSi Particular attention has been de- 
voted to atomic-scale devices, since they represent the 
limit towards one-dimensional electronics, and thus the 
transport properties of wire-like chains of atoms (espe- 
cially Au, Al and C), connected with metal electrodes 
have been widely investigatediiSiSSi^iS Despite their 
simple structures, atomic-sized chains display peculiar 
quantum properties due to their low dimensionality; in 
particular, the electronic properties are strongly affected 
by the nature of single chemical bonds and coordination 
numbers. 

The ongoing rapid advances in the measurements 
of electrical conductance in individual molecular- and 
atomic-sized devices require commensurate advances in 
the theoretical understanding of the detailed microscopic 



mechanisms that control charge mobility. Modeling of 
single nano-elements and coupled arrays of nano-devices 
is needed to provide interpretation and feedback to exper- 
imental measurements, to predict device characteristics, 
and to provide a basis for the functional progress of these 
new devices. 

In general, the electron transport properties of nanos- 
tructures can be simply described in the Landauer for- 
malism^ 3 *^ The Landauer relation connects the quan- 
tum conductance Q with the transmission function T(E) : 
its evaluation requires the knowledge of the electronic 
structure of the system under consideration and the in- 
clusion of scattering at contacts. This approach relies on 
the evaluation of lattice Green's functions of the system. 
Several approaches have been developed to calculate the 
quantum conductance in nanostructures, based on semi- 
empirical (tight- binding, Hiickel) models 25 26-2 ' 28 ; more 
recently, a variety of first-principles formulations have ap- 
pearedi 29 i 30 i 31 i 32 i 33 i 34 i 35 i 36 i 37 i 38 Ab-initio approaches have 
also been extensively used to characterize the electrical 
properties of nano- and bio-materialsf 3 ^* 4 ^ and to study 
the effects of microscopic structural relaxation and of 
electrode/conductor junctions^ 1 - 

In this paper, we present an original approach to the 
calculation of coherent transport properties of nanostruc- 
tures from first principles. Our methodology combines 
an accurate description of the electronic ground state 
provided by well-developed first-principles calculations 
based on plane-wave (PW) representations, with the Lan- 
dauer approach to describe transport properties of ex- 
tended systems^ 2 * 4 ^ The essential connection is provided 
by the use of the maximally-localized Wannier function 
representation— that allows to introduce naturally the 
ground-state electronic structure into the lattice Green's 
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FIG. 1: Geometry of a typical conductance calculation. 
Charge carriers can be transmitted throughout the contact 
region (grey) from the left to the right lead. Dotted lines sepa- 
rate the device in the three regions (left lead/conductor/right 
lead) that enter the conductance calculation. While the leads 
can be easily described with a primitive periodic unit cell, the 
conductor region has to be described with a supercell calcu- 
lation that includes - along the direction of electron trans- 
mission - both the conductor and the lead-conductor contact 
region (see Sec. IV). 



function approach that will be our basis for the evalua- 
tion of the Landauer quantum conductance. 

The paper is organized as follows. In Section [H] we 
describe the main features of the method. In Section ITTTI 
we will study the bulk electronic and conduction prop- 
erties of linear chains of Al and C atoms. Section IIVI 
deals with the conduction properties of a carbon nan- 
otube with a substitutional defect. This example is used 
to elaborate on the formulation of the two-terminal con- 
ductance problem in our approach. The paper ends with 
Section Ivl where our conclusions are presented. 



II. METHOD 
A. Electronic transport in extended systems 

Calculations of the quantum conductance are based 
on a recently developed efficient method for evaluat- 
ing quantum transport in extended systems^*^^ This 
method is applicable to any Hamiltonian that can be ex- 
panded within a localized-orbital basis and can be used 
as a general theoretical scheme for the computation and 
analysis of the electrical properties of nanostructures. 

Let us consider a system composed of a conductor con- 
nected to two semi-infinite leads, as in Fig.QJ The quan- 
tum description of the electronic conductance is a com- 
plex non-equilibrium problem. We begin the study of 
conduction properties focusing on the coherent electron 
transport. This approach leaves out non-equilibrium ef- 
fects due, e.g., to dissipative scattering or to an external 
bias. 

Quantum conductance (G) is the microscopic quantity 
that characterizes the transport properties of a conductor 
and may be calculated using the Landauer expression: 23 



Q(E) = ^T(E), 



where T(E) is the transmission function, i.e. the prob- 
ability that an electron injected at one end of the con- 
ductor with energy E will be transmitted to the other 
end. In the Green's function formalism the transmission 
function T is expressed as^ 



T = Tv(T L G c T R G c ), 



(2) 



where G^' a ^ are the retarded and advanced Green's func- 
tions of the conductor, and Ts^ m represent the coupling 
of the conductor to the leads. 

The effect of the semi-infinite leads on the conductor 
can be described by means of finite-dimension operators 
known as the self-energies The Green's function 

Gc and the coupling functions T^^y of Eq. [2] are ex- 
plicitly obtained from the self-energies as: 



Gc — (e — He — £l — S; 



{L,R} 



(3) 
(4) 

where He is the Hamiltonian matrix of the conductor, 
calculated with respect to a localized real-space basis set; 
in our approach all these operators will be represented 
in the basis set of the maximally-localized ground-state 
Wannier functions (WFs) for the system examined. This 
will allow a mapping for the first-principles Hamiltonian 
onto an exact band-structure tight-binding model, for 
which we then follow the detailed Landauer formulation 
of Ref. |H 

The key assumption is the description of the original 
system as a sequence of principal layers, 45 along which we 
calculate the transport properties and where the interac- 
tion between adjacent layers is accounted for by transfer 
matrices Tl,.r and Tl^r. Within this approach the self- 
energies due to the semi-infinite leads are 2 ^ 



£ x = H{ c {e - - {H^T L )-^H LC , 
= Hcr(c - H§) - H^ 1 T B y 1 H CR , 



(5) 



(1) 



where Hlc, Hcr are the coupling matrices for the 
conductor-lead assembly, and are the matrix ele- 

ments of the Hamiltonian for the infinite (bulk- like) leads. 

L R, 

For instance, H q describes the intralayer interactions 
and Hq{ R the interlayer coupling between two adjacent 
layers. If the orbitals are sufficiently localized, the resid- 
ual coupling with layers farther apart will be negligi- 
ble; conversely, the localization properties of the orbitals 
spanning the desired energy window determine the mini- 
mal thickness for a principal layer. The transfer matrices 
T l r and Ti, r are also computed from the Hamiltonian 
matrix elements via an iterative procedure^ 

The only required inputs then are the matrix elements 
of the Hamiltonian H mn expanded in a localized-orbital 
basis; by choosing the maximally-localized WFs repre- 
sentation, we provide essentially an exact mapping of the 
ground state onto a minimal basis4£ The accuracy of the 
results directly depends on having principal layers that 
do not couple beyond next-neighbors, i.e. on having a 
well-localized basis. 



B. Ab-initio electronic structure 

The starting point for our procedure is the first- 
principles calculation of the electronic structure of the 
nanostructure, eventually coupled to the leads. We adopt 
a standard electronic-structure method based on self- 
consistent total energy and force minimization, which al- 
lows to optimize simultaneously the atomic positions and 
the corresponding electronic wavefunctions. The elec- 
tronic structure is described within Density Functional 
Theory (DFT)4& The examples presented in the follow- 
ing Section are obtained in the Local Density Approxima- 
tions (LDA)4£ but more sophisticated corrections at the 
exchange-correlation functional {e.g. generalized gradi- 
ent approximation) can obviously be used. The electron- 
ion interaction is described via norm-conserving pseu- 
dopotentials 5 ^ in the form of Kleinman and Bylander. 51 
We use Periodic Boundary Conditions (PBC) along the 
three directions of direct space and the electronic wave- 
functions are expanded in a plane-wave basis set com- 
patible with the chosen PBC. While the translational in- 
variance makes the plane waves a very natural choice 
to describe the wavefunctions of a periodic system, the 
drawback is that they are truly delocalized. Brillouin 
Zone (BZ) summation are performed with homogeneous 
Monkhorst-Pack grids 52 of k-points in the the first Bril- 
louin Zone. 

This approach to electronic-structure calculations is 
widely used, provides a faithful description of the elec- 
tronic properties of the systems of interest, and in the 
present context it has been successfully applied to the 
investigation of both solid-state and molecular assem- 
bliesiS* 5 ^ The results of such calculations are, at each 
k-point, the Kohn-Sham energy eigenvalues and their 
corresponding eigenvectors (Bloch states), expanded in 
plane waves. It is worth noting that the present method- 
ology to compute quantum transport from first-principles 
will apply to any electronic-structure approach, since it 
can construct and employ orbitals that are maximally- 
localized, and that represent a minimal basis set, under 
the single assumption that eigenstates are in the Bloch 
form. Even if the basis set used is already localized (al- 
beit not minimal) the localization procedure will allow 
to recover exact results with smaller systems. In partic- 
ular, our procedure can be applied in combination with 
Car-Parrinello molecular dynamics simulations, opening 
the way to compute quantum conductance in large-scale 
systems, and as a function of temperature, completely 
from first-principles^ 5 - 



C. Maximally-localized WFs 

Bloch orbitals cannot be used directly to evaluate elec- 
tronic transport with the method outlined in Sec. Ill Al 
As we have pointed out, the quantum conductance is 
computed starting from the knowledge of the lattice 
Green's function, whose calculation relies on a local- 



ized orbital representation of the electronic states in real 
space. Bloch orbitals, that are intrinsically delocalized, 
have to be transformed into localized functions in order 
to construct the sparse, short-ranged matrix elements of 
the Hamiltonian. The core of our proposed methodology 
is to use maximally-localized WFs for the system con- 
sidered. These are the most natural choice for a set of 
localized orbitals that still span the same Hilbert space of 
the Hamiltonian eigenfunctions, and they allow to bridge 
plane-wave electronic structure and lattice Green's func- 
tion calculations in a coherent fashion. In the case of 
an isolated system the maximally-localized WFs become 
Boys localized orbitals^ therefore, our procedure is not 
tied to an extended-systems formulation, but can equally 
well represent isolated molecules. (In addition, the local- 
ization procedure is greatly simplified for the case of large 
unit cells, when T-sampling only is used—). 

A Wannier function w„R.(r), labeled by the Bravais 
lattice vector R, is usually defined via a unitary trans- 
formation of the Bloch functions V'nk(r) of the nth band: 



V 



(2tt) 2 



-ik R j3 



(6) 



BZ 



where V is the volume of the unit cell and the integration 
is performed over the entire Brillouin Zone. It is easy to 
show that the WFs defined as above form an orthonormal 
basis set, and that any two of them, for a given index n 
and different R and R', are just translational images of 
each other. Note that, as the Af WFs form a (continuous) 
linear combinations of Bloch functions with different en- 
ergies, they do not represent stationary states, but still 
span exactly the same original Hilbert space. The ab- 
initio eigenstates are well-defined, modulus an arbitrary 
k-dependent phase factor; thus, the definition above does 
not lead to a unique set of Wannier functions^, since the 
electronic structure problem is invariant for the transfor- 
mation tpnk ~* e^^i/Vik- Besides this freedom in the 
choice of phases <j) n (k) for the Bloch functions, there is a 
more comprehensive gauge freedom stemming from the 
fact that the many-body wavefunction is actually a Slater 
determinant: a unitary transformation between orbitals 
will not change the manifold, and will not change the to- 
tal energy and the charge density of the system. In all 
generality, starting with a set of N Bloch functions with 
periodic parts u n k, we can constructs infinite sets of J\f 
WFs displaying different spatial characteristics: 



V 



(2tt) 3 



BZ 



$>^™k(r) 



-ik-R j3 



d A k. (7) 



The unitary matrices include also the gauge freedom 
on phase factors afore mentioned 4» 

For our purposes, we need to transform the Bloch 
eigenstates in WFs with the narrowest spatial distribu- 
tion. We construct maximally-localized WFs using the 
algorithm proposed by Marzari and Vanderbilt44 We 
define a Spread Operator (f2) as the sum of the second 



moments of the Wannier functions corresponding to one 
choice of translational lattice vector: 

fi = ^[(w n o\r 2 \w n0 ) - (w„o|r|w Il0 ) 2 ], (8) 

n 

where the sum is over the group of bands which spans 
the Hilbert space. The value of the spread depends 
on the choice of unitary matrices U^; thus it is possible 
to evolve any arbitrary set of until the minimum 
condition 



is satisfied. At the minimum, we obtain the matrices 
(Jj( k )) ML that transform the first-principles V^if ( r ) into 
the maximally-localized WFs (r) : 

(10) 

A useful feature of the method is that the only ingre- 
dients needed to calculate the spread functional Q and to 
evolve the unitary matrices are the overlap matrix 
Mmn 3 ^ between the periodic part of the Bloch states at 
neighboring k-points: 

M mn ] = (w m ,kK,k+b>, (11) 

where b is the vector that links neighboring k-points in 
the discretized BZ integrals^ 

It is important to notice that whenever a Born-von 
Karman discretization of the Brillouin Zone is intro- 
duced, even the above-mentioned WFs are not truly 
localized, but will be periodic in real-space, with a 
superperiodicity determined by the BZ discretization. 
The truly isolated limit is recovered only in the case of 
continuous BZ integrations. This is easily seen remem- 
bering that V'nk(r) = u„k(r)e ik r , and w„k(r) has the 
periodicity of the direct lattice; thus the phase factors 
e ik-r determine the superperiodicity of the ^„k them- 
selves. In the standard language of electronic-structure 
calculations, if the ip n k have k's that are restricted to a 
uniform Monkhorst-Pack mesh, they will all be periodic 
with a wavelength inversely proportional to the spacing 
of the mesh; this periodicity is consequently inherited by 
the WFs. For Af k-points along a direction of the BZ, 
the WFs will repeat along the corresponding direction 
every Af cells; therefore a mesh of k-points needs to 
be dense enough to assure that adjacent replicas of the 
WFs do not overlap. 

The method described above works properly in the 
case of isolated groups of bandsi^S On the other hand 
to study quantum conductance in extended systems we 
often need to compute WFs for a subset of energy bands 
that are entangled or mixed with other bands. Most of- 
ten we are interested in the states that lie in the vicinity 



of the Fermi level of a conductor in a restricted energy 
range. Since the unitary transformations mix energy 
bands at each k-point, any arbitrary choice of states in- 
side a prescribed window will affect the localization prop- 
erties of WFs unless energy gaps effectively separate the 
manifold of interest from higher and lower bands. This 
problem has been solved by Souza, Marzari, and Van- 
derbilt, introducing an additional disentanglement pro- 
cedural that automatically extracts the best possible 
manifold of a given dimension from the states falling in 
a predefined energy window. This is the generalization 
to entangled or metallic cases of the maximally-localized 
WF formulation. The procedure relies on minimizing the 
subspace dispersion across the Brillouin Zone, and effec- 
tively extracts the bands of interest from the overall band 
structure. In practice, first we select a desired number 
of bands in an energy window; then we determine the 
optimally-connected subspace that can be extracted from 
that band structure; and finally we proceed with a stan- 
dard localization procedure inside the selected subspace, 
using the same kind of spread functional f2 and of unitary 
matrices Umn- The resulting orbitals have the same good 
localization properties, and allow to apply our formal- 
ism to arbitrary systems, independently of the insulating 
or metallic nature of the band manifold. It should be 
stressed that the WFs obtained in the later case are not 
the WFs of the occupied subspace (that would exhibit 
poor localization properties), but are those of a well con- 
nected, continuous subspace that in general will contain 
both occupied and unoccupied Bloch functions. 

In order to calculate the conductance according to the 
prescriptions outlined in Sec. Ill Al we need as an input 
the matrix elements of the Hamiltonian calculated on a 
localized basis: in our case, it is the minimal basis of 
the maximally-localized WFs. The advantages of this 
choice are twofold: firstly, besides being a minimal basis, 
the WFs span exactly the Hilbert space of an insulator 
and, with arbitrary accuracy, of an entangled metallic 
system. Secondly, their localization assures the choice of 
the system with the smallest number of atomic layers. 
The Hamiltonian matrices (H^, He, Hlc, Hqr) can 
be formally obtained from the on site (-ffoo) and coupling 
(Hqi) matrices between principal layers. In our formal- 
ism, and assuming a BZ sampling fine enough to elim- 
inate the interaction with the periodic images, we can 
simply compute these matrices from the unitary matrix 
£A k J obtained in the localization procedure^ By defini- 
tion of energy eigenvalues (e m k), the Hamiltonian matrix 
H mn (k) = eink^m,?!) is diagonal in the basis of the Bloch 
eigenstates. We can calculate the Hamiltonian matrix in 
the rotated basis, 

H^ rot \k) = (UWyH(k)uW. (12) 

Next we Fourier transform ff( rot )(k) into a set of Nk p 
Bravais lattice vectors R within a Wigner-Seitz supercell 



centered around R=0 : 

(13) 

where N^ p derives from the folding of the uniform mesh 
of k-points in the BZ. The term with R=0 provides the 
on site matrix Hqq = (w m o\H\w n o), an d the term R=l 
provides the coupling matrix Hqi = {w m o\H\w n i): These 
are the only ingredients required for the evaluation of the 
quantum conductance. 



III. BULK-LIKE CONDUCTANCE 

As a first application of our method, we consider a case 
in which leads and conductor (as sketched in Fig. ^) are 
made of the same material, and we compute the conduc- 
tance of the ideal and infinite nanostructure (bulk-like 
conductance). In this case, it is not necessary to distin- 
guish between conductor and lead terms and the single 
layer Hon and the coupling Hqi matrices are the only 
necessary input. 

We will focus on one-dimensional (ID) linear chains 
of atoms. The systems that have been studied most are 
chains of Au, 20 ! 21 ! 22 ! 26 . 63 ! 64 ! 65 , A&S1 and 0,68,89,7 71 ^ 
the following, we will discuss results for Al (Sec. IIII A|) 
and C (Sec. IIIIB|) chains. 



Aluminum chain 



— Al— Al — 




1 2 3 
G(2e 2 /h) 



FIG. 2: Linear Al-chain. Isosurface of (a) a a and (b) two n 
WFs. (c) Comparison between original first-principles (grey 
dots) and interpolated (black lines) band structure. Dotted 
lines represent the Fermi Level. Note that some bands (i.e. 
at Fermi Level) are double degenerate, (d) Quantum con- 
ductance plots. The solid line is relative to the whole energy 
range (E € [—7,6] eV); the shaded area identifies the quan- 
tum conductance calculated in the narrow energy window (E 
6 [—4.5, 1.5] eV) bounded by the dashed lines. 



An ideal and infinite Al chain is simulated using peri- 
odic boundary conditions and a unit cell containing two 
aluminum atoms. A large vacuum region (~ 10 A) in 
the direction perpendicular to the chain prevents the in- 
teraction with adjacent replicas. A (12x1x1) grid 
of k-points and 18 Ry energy cut-off for the wavefunc- 
tion expansion assure the convergence of the electronic 
structure of the system. The optimized Al-Al distance 
(d = 2.42 A) is in very good agreement with previous 
DFT investigations&Si 

Following the procedure described above, we calcu- 
lated the electronic structure and quantum conductance 
of this system. To construct the WFs we selected an 
energy window with E € [—7, 6] eV around the Fermi 
level (taken as the reference zero). This energy win- 
dow contains all the occupied bands and the first empty 
states. We chose to extract an 8-dimensional manifold 
from this energy window: After the disentanglement and 
localization procedure, we obtain eight WFs which span 
the 8-dimensional Hilbert manifold and represent an or- 
thonormal minimal basis for it. Thus, the calculation 
of quantum conductance involves operations with very 
small (8 x 8) matrices, with a negligible computational 
effort, exactly comparable to a tight-binding calculation 
(TB) with two sites and four orbital per site. However, 
our results provide more information on the electronic 



structure than the TB approach. The calculated WFs 
(Fig. |2K>b) are well characterized and are consistent with 
the estimated chemical bonds present in the system (cr 
and 7T orbitals). The a states are centered in the middle 
of the Al-Al bond, while the tt states are localized around 
single atoms. 

As a test of the accuracy of the WF transformation, we 
have computed back the band structure of the Al-chain, 
starting from the Wannier-function Hamiltonian in real 
space. This Hamiltonian can now be Fourier transformed 
in reciprocal space HmiP (R) (Eq. I13fl for any arbitrary 
k-point 

HkfW = ^e+ ik '- R ^)(R); (14) 

R 

the resulting Hamiltonian matrices can then be diagonal- 
ized to find energy eigenvalues. Comparing the original 
PW (grey dots) with the interpolated (black lines) band 
structure (Fig. |2t ) we see an excellent agreement. This 
is an expected but important validation, since it proves 
that the intermediate transformations do not affect the 
accuracy of the first-principles PW calculations. All the 
information on the electronic structure of the system is 
transferred to the matrix elements of the Hamiltonian 
expressed in the localized WFs basis. 



A closer inspection to Fig. |2t shows some unmatched 
electronic states at energies lower than the highest inter- 
polated band (in this case, in the vicinity of the T point). 
This is the effect of the band-space minimization, 61 which 
singles out the best-connected manifold from an entan- 
gled group of states. The upper edge of the outer win- 
dow (Fig. intersects states having comparable energy 
but different symmetries, relevant to higher parts of the 
spectrum. The inclusion of these contributions (i.e., the 
states around T at ~ 4-5 eV) would affect significantly 
the localization properties of the WF basis. The mini- 
mization of the dispersion for the extracted manifold (the 
disentanglement of the bands) is thus an essential step in 
the WFs calculation. 

The disentanglement procedure can be used to probe 
different energy windows; this allows to single out the 
most relevant bands, linking conductance properties to 
the nature of the chemical bonds. As an example, if 
we restrict the energy window to a few eV around the 
Fermi Level (dashed lines in Figure |2t,d), we can de- 
scribe quantum conductance using the only smaller set 
of bands (three) included in the narrow window. In this 
restricted range, the new conductance spectra (shaded 
area in Fig. EJl) is indistinguishable from the original one 
(black line). 

The linear chain of Al atoms displays metallic behav- 
ior (Fig. d), in agreement with previous DFT calcula- 
tions^ Due to the reduced coordination number of the 
Al atoms in the chain (n c hain = 2) compared with the 
FCC-bulk phase (ribulk = 12), this metallic character was 
not obvious a priori. It is important to note that, in 
general, metallic systems are not well represented in a 
WF framework. So far, only transition metals have been 
the subject of WFs studies, due to the localized charac- 
ter of their d-orbitalsj22 and only recently localized wave 
functions in reciprocal space have have been proposed 
for simple metals (Na and Al)iS However, in the dis- 
entanglement procedure we are not required to restrict 
ourselves to the occupied subspace, but we can mix filled 
and empty states, allowing us to extract well-connected 
manifolds that have the same localization properties of 
the manifolds for insulators and semiconductors. 

To better understand the conductance properties of 
Al chains, we have also calculated the eigenvectors of 
the transmission function T(E), generally known as 
eigenchannels*^ The eigenchannels completely charac- 
terize transmission, and at each energy, describe the sin- 
gle modes of the electronic transport. Our results show 
that the eigenchannels at the Fermi energy are barely the 
linear combination of the 7r-like WFs of Figure^- The 
two quanta of conductance (2e 2 /h) at the Fermi level 
(Fig. 0i) correspond to two degenerate it states, which 
constitute the channels for charge mobility. The metallic 
behavior is in qualitative agreement with the geometrical 
properties of WFs. As mentioned before, both a- and n- 
WFs are well localized. While the a states (Fig. are 
centered in the middle of the Al-Al bond, in a bonding 
configuration, the 7r orbitals (Fig.^Jj), responsible for the 



chain metallicity, are centered on the single atoms. 

Finally, in the absence of external leads, there is a 
one-to-one correspondence between the quantum conduc- 
tance spectrum (Fig. 2d) and the band structure: at a 
given value of E, the quantum conductance (Eq. ^| is 
a constant proportional to the number of transmitting 
channels available for charge mobility, which are equal (in 
a periodic system) to the number of bands at the same 
energy. The perfect agreement between band structure 
and quantum conductance represents a further valida- 
tion for the ability of our method to calculate transport 
properties. 



B. Carbon chains 

As a second application, we have studied two differ- 
ent species of carbon chains. Nanodevices where C- 
chains act as conductors bridging metal electrodes are 
not only ideal prototypes for studying conduction in 
reduced-dimensionality systems, but are also fundamen- 
tal constituents of low-pressure carbon assemblies, such 
as those found in end-capped molecules or in the inter- 
stellar medium. Theoretical models proposed so far have 
dealt with wires of equidistant C atoms trapped between 
metallic leads of Au or Ali^LS In those cases, the con- 
duction properties of the system are strongly dependent 
on the number (odd or even) of atoms in the chain. 

Here, we focus on the effects of structural relaxation on 
the electronic and transport properties of infinite carbon 
chains, known as carbyne. The name carbyno&l&H de- 
notes an allotrope based on a linear chain of sp-hybridized 
carbon atoms: isomeric polyethynylene diylidene (poly- 
cumulene or cumulene) or polyethynylene (polyyne) . The 
cumulene form is characterized by an equidistant ar- 
rangement of C-atoms with double sp-bonds (= C = C = 
)„, while the polyyne form is a dimerized linear chain with 
alternating single-triple bonds (— C = C—) n . The exper- 
imental evidence for carbyne chains is controversial and 
its properties not completely known. 68 1 76 We studied the 
effects of the two allotropes (cumulene vs. polyyne) to 
the electronic and conduction properties of carbyne. 

We used four C-atoms in a periodically repeated cell, 
and an (8 x 1 x 1) grid of k-points for BZ summation. 
The electronic wavefunctions are expanded in a plane- 
wave basis set up to 40 Ry. We first optimized the lat- 
tice constant of the cumulene structure, and then, in 
the same unit cell, we relaxed the carbon-carbon dis- 
tances in the polyyne phase. In the cumulene form the 
C-atoms are separated by d cumu i ene = 1.37 A, while 
the polyyne form dimerizes with C-atoms separated by 
^■single = 1-51 A and dtripie = 1-22 A, in agreement with 
previous theoretical calculations. — Incidentally, polyyne 
is energetically more stable than cumulene by 1.2 eV per 
unit (C-C). 

Fig. |21 and Fig. 0] show our results. Cumulene 
(Fig. b) is characterized by symmetric sp-bonds, uni- 
formly distributed along the chain, a states are localized 
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FIG. 3: Cumulene. Isosurface of (a) two a and (b) two n WFs. 
First-principles band structure (c), and calculated quantum 
conductance (d) in the selected energy window. Dotted lines 
represent the Fermi Level of the system. 
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FIG. 4: Polyyne. Isosurface of (a) two a and (b) two n WFs. 
First-principles band structure (c), and calculated quantum 
conductance (d) in the selected energy window. Dotted lines 
represent the Fermi Level of the system. 



in the middle of C-C bonds while it states are centered 
around single C-atoms. In polyyne (0Ji,b), a orbitals are 
localized both on single C—C and on triple C = C bonds, 
with a a state in the middle of each bond. The tt orbitals 
are localized on the C = C bonds: there are two of these 
7r orbitals in the middle of each triple bond, related by a 
90° rotation around the axis. 

The electronic structures and quantum conductances 
of cumulene and polyyne are shown in Fig. |2t,d and 
Fig- Etjd respectively. The symmetric chain displays 
metallic behavior, in agreement with previous theoret- 
ical calculations^! Polyyne is instead semiconducting: 
the relaxation of the carbon-carbon distances induces a 
Peierls-type distortion, which stabilizes the structure and 
opens energy gaps at the edges of the Brillouin Zone. 
The metallicity of cumulene is an effect related to the 
homogeneous distribution of the atoms, and not to the 
dimensionality of the chain: polyyne, which has the same 
dimensionality of cumulene, is not a metal. As men- 
tioned in the previous section, the electrical (metallic or 
semiconducting) behavior is tightly reflected in the geo- 
metrical properties of WFs. The eigenchannels of both 
systems, near the Fermi energy, are made by linear com- 
bination of the 7r-like WFs of Figures 1314b . As was the 
case for the aluminum chain, the n orbitals of cumulene 
are located on the atoms, and the system is metallic. 
Contrary-wise, in polyyne the it states are centered in 
the middle of the triple C-bonds, and the system is semi- 
conducting. 



IV. TWO TERMINAL CONDUCTANCE 

As a prototypical example of a two terminal device 
we have computed the quantum conductance for a small 
zigzag (5,0) carbon nanotube in the presence of an iso- 
lated substitutional silicon defect. Topological and sub- 
stitutional defects in carbon nanotubes have been widely 
studied^2*22iSi2S and our results can be directly com- 
pared with available theoretical data. 

Using Fig. ^ a s reference, we choose in the present 
case as the conductor region a finite segment of nanotube 
which contains the defect, while the leads are modeled by 
two semi-infinite nanotubes. To compute the two termi- 
nal conductance we need, in principle, three sets of cal- 
culations (two only if the leads are of the same material) : 
bulk calculations for the two infinite leads and a super- 
cell calculation for the conductor and the contacts (see 
Fig. ^) . The supercell needs to be chosen large enough 
so that the influence of the conductor wavefunctions on 
the leads becomes negligible at the boundaries, assuring 
seamless matching between the lead WFs determined in 
the supercell and in the bulk calculations. The match- 
ing condition can also be expressed by saying that the 
on-site and hopping integrals (in the tight-binding lan- 
guage) of the individual WFs have to be the identical on 
both sides of the interface boundary. The matrices H ' 
and Hq '^ for the infinite leads are obtained from the 
bulk calculations (in this case, the ideal (5,0) nanotube), 
while the other coupling matrices are derived from the 
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FIG. 5: a) Isosurface of two a WFs in a (5,0) carbon nanotube 
in presence of a substitutional Si defect (black atom). The 
polarization of the a states is due to the effect of the Si defect, 
b) Quantum conductance plots for the (5,0) nanotube with 
and without the Si defect (top panel); same nanotube with a 
periodic arrangement of Si defects in a bulk-like configuration 
(bottom panel). Vertical dotted lines represent the Fermi 
Level of the systems. 

supercell calculation. We stress the need to include a 
sufficiently large portion of the contacts in the supercell 
calculation so that, far from the interface and into the 
leads, the system recovers its bulk properties. The char- 
acteristic length determining the amount of conductor 
included is not the decay length of the density matrix of 
the metal (long-ranged, and algebraic), but the charac- 
teristic length of the WFs disentangled from the metal, 
and whose localization properties on these well-connected 
manifold are comparable to those of a semiconductor or 
an insulator. The properties of uniqueness and localiza- 
tion of the WFs ensure effortless continuity across the 
interface into the bulk leads. Note that a similar defini- 
tion of the geometry of the simulation cell is needed also 
when using other real-space methods to compute quan- 
tum conductancei^Si^ Since our supercell contains both 
the conductor and a portion of the leads large enough 
to take into account the presence of the contacts, the 
Wannier transformation produces a set of WFs cover- 
ing the whole coupled region. However, since the WFs 
are strongly localized, it is straightforward to distinguish 
those centered on the conductor and those on the leads, 
and to see the matrix elements of the Hamiltonian seam- 



less turn away from the conductor into bulk-like lead 
terms. This approach brings another advantage: since 
the conductor and the conductor-lead interfaces are sim- 
ulated in the same supercell, we have a consistent (and 
fully relaxed) microscopic picture of the contacts. 

The (5,0) carbon nanotube (see Fig. has been 

simulated in a (4.25 x 10.0 x 10.0) A 3 supercell, with 8 k- 
points and a 40 Ry cutoff for the wavefunction expansion. 
Due to the highly localized nature of the WFs, only four 
atomic layers are needed to reproduce the bulk-like be- 
havior of the nanotube far away from the Si defect. The 
presence of the defect does not induce significant struc- 
tural relaxations, but affects the electronic and conduc- 
tion properties of the tube. Figure shows the isosur- 
faces of two calculated WFs with a character. Far away 
from the defect, WFs are symmetric and centered in the 
middle of the C-C bond. Closer to the defect, the differ- 
ent electronegativity of silicon and carbon atoms induces 
a polarization of the Si-C bond that modifies the con- 
duction properties of the system. Figure |3Jd (top panel) 
displays the quantum conductance of the nanotube with 
respectively in the presence of a defect (C(5,0)/Si, solid 
line) or in the ideal clean case (C(5,0), dashed line). The 
conductance of the ideal (5,0) nanotube shows the typ- 
ical step- like shape and a metallic behavior, which is 
the result of the high curvature of such a small radius 
nanotubei 81 ! 82 Once a single Si defect is introduced, the 
system maintains its metallic character but the overall 
spectrum changes drastically. We observe a general re- 
duction of the conductance along the whole energy range, 
and the distortion of the step-like shape of the pure nan- 
otube. The appearance of dips, corresponding to the 
discontinuities in the original step function, is a char- 
acteristic feature of nanotubes with defects (SLSiSl where 
the backscattering of electrons reduces the quantum con- 
ductance. 

As a consequence of the external leads the quantum 
conductance is not directly related to the band structure 
of the supercell. We show this in Fig.fSjD (bottom panel), 
where the bulk conductance of the C(5,0)/Si system is 
presented. The general trend of (5(E) is different from 
both curves in the upper panel. The whole C(5,0)/Si 
system is periodically repeated, and the conductance as- 
sumes the typical step-like behavior. However, with re- 
spect to the pure nanotube, we observe both an overall 
increse of the quantum conductance and the presence 
of novel features, such as the marked peak due to the 
silicon states, just below the Fermi level. The periodic 
distribution of the Si atoms leads to the formation of ad- 
ditional channels available for the charge transport along 
the tube. On the other hand, these features disappear 
in the case of the single defect, where the breaking of 
the traslational symmetry does not allow the formation 
of delocalizcd orbitals. In conclusion, while the doping 
with a regular pattern of Si-atoms increases the conduc- 
tance, the scattering dephasing of a single defect reduces 
the global transport properties of the small carbon nan- 
otube. 



V. CONCLUSIONS 

We have presented a novel approach to calculate quan- 
tum conductance in extended systems in the coherent 
transport regime. Our methodology combines the accu- 
rate self-consistent minimization of the ground-state elec- 
tronic structure via first-principles calculations, the de- 
termination of the maximally-localized WFs correspond- 
ing to manifold of bands spanning the energy range rele- 
vant for conduction, and the calculation of the quantum 
conductance using a real-space Green's function formal- 
ism based on the Landauer approach. This procedure 
opens the way to selectively describe the quantum con- 
ductance in terms of the relevant one-electron states that 
contribute directly to the transport process. It links the 
description of electronic conductance to the intrinsic na- 
ture of the chemical bond, and gives new insights onto 
the essential mechanisms that govern the electron flow 
at the nanoscale. Moreover, it is computationally very 
efficient and can be straightforwardly implemented as a 
post-processing step in any standard electronic structure 
calculation 83 leading to a first-principles, highly accurate 
computation of electron transport properties. 



As a first illustration of the potential of this method- 
ology, we have studied quantum conductance in linear 
chains of aluminum and carbon atoms, and in defective 
carbon nanotubcs. In all cases, we have underlined the ef- 
fects of the reduced coordination and of the atomic relax- 
ation on the transport properties and we have established 
a clear relationship between the electrical characteristics 
and the chemical bonds in the system. 
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